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SYSTEMS FOR GENERATING AND USING A LOOKUP TABLE WITH 
PROCESS FACILITY CONTROL SYSTEMS AND MODELS OF THE SAME, 
AND METHODS OF OPERATING SUCH SYSTEMS 

COPYRIGHT NOTTPF 

A portion of the disclosure of this patent document (software listings in 
Appendices A and B) contains material that is subject to copyright protection. The 
copyright owner has no objection to the facsimile reproduction by anyone of this patent 
document or the patent disclosure, as it appears in the United States Patent and 
Trademark Office patent file or records, but otherwise reserves all copyright rights and 
protection whatsoever. 



CROSS-REFER FNCE TO RELATED PATENT nnr-rr^Ni TS 

The present invention is related to that disclosed in (i) United States Patent No. 
5,351,184 entided "Method of Multivariable Predictive Control Utilizing 
Range Control;" (ii) United States Patent No. 5,561,599 entitled "Method of 
Incorporating Independent Feedforward Control in a Multivariable 
PREDICTIVE CONTROLLER;" (iii) United States Patent No. 5,574,638 entitled "METHOD 
OF OPTIMAL SCALING OF VARIABLES IN A MULTIVARIABLE PREDICTIVE CONTROLLER 

Utilizing Range Control;" (iv) United States Patent No. 5,572,420 entitled 

"METHOD OF OPTIMAL CONTROLLER DESIGN OF MULTIVARIABLE PREDICTIVE CONTROL 

Utilizing Range Control" (the '"420 Patent"); ( v ) United States Patent No. 
5,758,047 entitled "Method of Process Controller Optimization in a 
Multivariable Predictive Controller;" (vi) United States Patent Application 
Serial No. 08/490,499, filed on June 14, 1995, entitled "Method of Process Controller 
Optimization in a Multivariable Predictive Controller;" (vii) United States Patent 
Application Serial No. 08X850,288 entitled "Systems And Methods For Globally 
Optimizing a Process Facility;" (viii) United States Patent Application Serial No. 
08X851,590 entitled "Systems for Generating and Using a Lookup Table with 
Process Facility Control Systems and Models of the Same, and Methods of 
Operating Such Systems;" and (ix) United States Patent Application Serial N 7 
09X137,358 entitled "Controllers That Determine Optimal Tuning Parameters 
For Use in Process Control Systems And Methods of Operating The Same;" and 
(x) United States Patent Application Serial No. (Attorney Docket No. 120 25206), 
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entitled "Process Facility Control Systems Using an Efficient Prediction Form 
And Methods of Operating The Same" (which application is filed concurrently 
herewith), all of which are commonly assigned to the assignee of the present invention. 
The disclosures of these related patent applications are incorporated herein by reference 
5 for all purposes as if fully set forth herein 

TECHNICAL FIELD OF TP F. INVENTION 
The present invention is directed, in general, to control systems for process 
facilities and, more specifically, to systems for generating and using lookup tables with 
10 process facility control systems and models of the same, and methods of operating such 
systems, all for use to optimize process facilities. 

BACKGROUND OF THE I NVENTION 

Presently, process facilities (e.g., a manufacturing plant, a mineral or crude oil 
15 refinery, etc.) are managed using distributed control systems. Contemporary control 
systems include numerous modules tailored to control or monitor various associated 
processes of the facility. Conventional means link these modules together to produce 
the distributed nature of the control system. This affords increased performance and a 
capability to expand or reduce the control system to satisfy changing facility needs. 

Process facility management providers, such as Honeywell, Inc., develop 
control systems that can be tailored to satisfy wide ranges of process requirements (e.g., 
global, local or otherwise) and facility types (e.g., manufacturing, refining, etc.). A 
primary objective of such providers is to centralize control of as many processes as 
possible to improve an overall efficiency of the facility. Each process, or group of 
associated processes, has certain input (e.g., flow, feed, power, etc.) and output (e.g., 
temperature, pressure, etc.) characteristics associated with it. 

In recent years, model predictive control ("MPC") techniques have been used to 
optimize certain processes as a function of such characteristics. One technique uses 
algorithmic representations to estimate characteristic values (represented as parameters, 
variables, etc.) associated with them that can be used to better control such processes^ 
In recent years, physical, economic and other factors have been incorporated into 
control systems for these associated processes. Examples of such techniques are 
described in United States Patent No. 5,351,184 entitled "Method of Multivariable 
Predictive Control Utilizing Range Control;" United States Patent No. 
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5,561,599 entitled "Method of Incorporating Independent Feedforward 
Control in a Multi variable Predicts Controller;" United States Patent No 
5,574,638 entitled "Method of Optimal Scaling of Variables in a Multi variable 
PRED,CTIVE Controller Utilizing Range Control;" United States Patent No 
5,572,420 entitled "Method of Optimal Controller Design of Multtvariable 
Predictive Control Utilizing Range Control" (the "<420 Patent"); United States 
Patent Application Serial No. 08\850,288 entitled "Systems And Methods For 
Globally Optimizing a Process Facility;" United States Patent Application Serial 
No. 08\851,590 entitled "Systems and Methods U SI ng Bridge Models to 
Globally Optimize a Process Facility;" and United States Patent Application Serial 
No. 09U37.358 entitled "Controllers That Determ^e Optimal Tuning 

PARAMETERS FOR USE IN PROCESS CONTROL SYSTEMS AND METHODS OF OPERATING 

THE SAME," all of which are commonly owned by the assignee of the present invention 

and incorporated herein above by reference for all purposes. 

Generally speaking, one problem is that conventional efforts, when applied to 

specie processes, tend to be non-cooperative (e.g., non-global, non-facility wide etc ) 
and may, and all too often do, detrimentally impact the efficiency of the process facility 
as a whole. For instance, many MFC techniques control process variables to 
predetermined set points. Oftentimes the set points are a best estimate of a value of the 
set pomt or set points. When a process is being controlled to a set point, the controller 
may not be able to achieve the best control performances, especially under 
process/model mismatch. 

To further enhance the overall performance of a control system, it is desirable to 
design a controller that deals explicitly with plant or model uncertainty The <420 
Patent, for example, teaches methods of designing a controller utilizing range control 
The controller is designed to control a "worst case" process. An optimal controller for 
the process is achieved and, if the actual process is not a "worst case process " the 
performance of the controller is better than anticipated. 

There are a number of well known PID "tuning" formulas, or techniques and 
the most common, or basic, PID algorithm includes three known user specified tuning 
parameters (K, 2) whose values determine how the controller will behave ThesT " 
parameters are determined either by trial and error or through approaches that require 
knowledge of the process. Although many of these approaches, which are commonly 
algorithms, have provided improved control, PID controller performance tuned by such 
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algorithms usually degrades as process conditions change, requiring a process engineer, 
or operator, to monitor controller performance. If controller performance deteriorates, 
the process engineer is required to "re-tune" the controller. 

Controller performance deteriorates for many reasons, although the most 
common cause is changing dynamics of the process. Since PID controller performance 
has been related to the accuracy of the process model chosen, a need exists for PID 
controllers that allows for such uncertainty by accounting for changing system 
dynamics. Further, the requirement for ever-higher performance control systems 
demands that system hardware maximize software performance. Conventional control 
system architectures are made up of three primary components: (i) a processor, (ii) a 
system memory and (iii) one or more input/output devices. The processor controls the 
system memory and the input/output ("I/O") devices. The system memory stores not 
only data, but'also instructions that the processor is capable of retrieving and executing 
to cause the control system to perform one or more desired functions. The I/O devices 
are operative to interact with an operator through a graphical user interface, and with the 
facility as a whole through a network portal device and a process interface. 

Over the years, the quest for ever-increasing process control system speeds has 
followed different directions. One approach to improve control system performance is 
to increase the rate of the clock that drives the system hardware. As the clock rate 
increases, however, the system hardware's power consumption and temperature also 
increase. Increased power consumption is expensive and high circuit temperatures may 
damage the process control system. Further, system hardware clock rate may not 
increase beyond a threshold physical speed at which signals may be processed. More 
simply stated, there is a practical maximum to the clock rate that is acceptable to 

conventional system hardware. 

An alternate approach to improve process control system performance is to 
increase the number of instructions executed per clock cycle by the system processor 
("processor throughput"). One technique for increasing processor throughput calls for 
the processor to be divided into separate processing stages. Instructions are processed 
in an "assembly line" fashion in the processing stages. Each processing stage ^ 
optimized to perform a particular processing function, thereby causing the processor as 
a whole to become faster. There is again a practical maximum to the clock rate that is 
acceptable to conventional system hardware. 

Since there are discernable physical limitations to which conventional system 
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hardware may be utilized, a need exists broadly for an approach that decreases the 
number of instructions required to preform the functions of the process control system. 
A need exists for such an approach that accounts for process uncertainty by accounting 
for changing process dynamics. 

SUMMA RY OF THF, INVENTION 
To address the above-discussed deficiencies of the prior art, it is a primary 
object of the present invention to provide systems and methods of operating such 
systems for populating and using lookup tables with process facility control systems, as 
well as models of the same. In accordance with an exemplary embodiment below- 
discussed, the principles of the present invention may be used to define and populate a 
lookup table in response to the needs of a global controller. The lookup table is 
populated with a range of possible values of at least one measurable characteristic 
associated with one or more processes of the process facility and in accordance with a 
model of at least a portion of the same. 

Rather than calculate and re-calculate certain characteristics associated with a 
process or process model, which would consume significant system resources, the 
present invention introduces a data structure capable of maintaining a range of possible 
values of one or more of such certain characteristics. Use of the lookup table in lieu of 
execution and re-execution of the instructions for performing characteristic calculations 
decreases the number of instructions required to preform the functions of the process 
control system. The lookup table, once suitably populated, accounts for process 
uncertainty by maintaining the range of possible values, thereby accounting for 
changing process dynamics. 

An exemplary computer system for use with a process facility that is capable of 
populating a data structure in accordance with the principles of the present invention 
includes both a memory and a processor. The memory is capable of maintaining (i) the 
data structure, which has a plurality of accessible fields, and (ii) a model of at least a 
portion of at least one process of a plurality of associated processes of the process 
facility. The model may advantageously include a mathematical representation of at 
least a portion of the at least one process, defining certain relationships among inputT 
and outputs of the at least one process. The processor is capable of populating ones of 
the plurality of accessible fields of the data structure using the model iteratively with a 
range of possible values of the at least one measurable characteristic. The computer 
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system is capable of using the range of possible values of the at least one measurable 
characteristic to predict an unforced response associated with the at least one process. 

In accordance with an important aspect hereof, the data structure may be 
populated and maintained on-line (e.g., at a controller, distributed through a process 
control system, etc.), off-line (e.g., standalone computer, computer network, etc.), or 
through some suitable combination of the same. Likewise, the data structure may 
remain static upon population, be dynamic, or be modifiable, at least in part. 

Those skilled in the art will understand that "controllers" may be implemented in 
hardware, software, or firmware, or some suitable combination of the same, and, in 
general, that the use of computing systems in control systems for process facilities is 
known. The phrase "associated with" and derivatives thereof, as used herein, may mean 
to include, be included within, interconnect with, contain, be contained within, connect 
to or with, couple to or with, be communicable with, cooperate with, interleave, be a 
property of, be bound to or with, have, have a property of, or the like; the term 
"include" and derivatives thereof, as used herein, are defined broadly, meaning 
inclusion without limitation; and the term "or " as used herein, means and/or. 

The foregoing has outlined rather broadly the features and technical advantages 
of the present invention so that those skilled in the art may better understand the 
detailed description of the invention that follows. Additional features and advantages of 
the invention will be described hereinafter that form the subject of the claims of the 
invention. Those skilled in the art should appreciate that they may readily use the 
conception and the specific embodiment disclosed as a basis for modifying or designing 
other structures for carrying out the same purposes of the present invention. Those 
skilled in the art should also realize that such equivalent constructions do not depart 
from the spirit and scope of the invention in its broadest form. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present invention, and the advantages 
thereof, reference is now made to the following descriptions taken in conjunction with 
the accompanying drawings, wherein like numbers designate like objects, and in which^ 

FIGURE la illustrates a simple block diagram of an exemplary process facility 
with which the present invention may be used; 

FIGURE lb illustrates a detailed block diagram of one of the exemplary local 
controllers introduced in FIGURE la; 
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FIGURE 2 illustrates a flow diagram of an exemplary method for populating a 
data structure in accordance with the principles of the present invention; and 

FIGURE 3 illustrates an exemplary two-dimensional graphical representation of 
MV and PV curves in accordance with the principles of the present invention. 

DETAILED DESCRIPTION OF THE DRAWINGS 

In accordance with the above-given summary, computer systems, and methods 
of operating the same, are introduced herein for populating and using lookup tables with 
process facility control systems, as well as models of the same. Before undertaking a 
detailed description of an advantageous embodiment of the present invention, and 
discussing the various benefits and aspects of the same, it is useful to understand 
conceptually the operation and control structure of an exemplary process facility. 

Initial reference is therefore made to FIGURE la, wherein a simple block 
diagram of such a process facility (generally designated 100) is illustrated. Exemplary 
process facility 100 is operative to process raw materials, and includes a control center 
105, six associated processes 110a to HOf that are arranged into three stages and a 
control system (generally designated 1 15). The term "include," as well as derivatives 
thereof, as used throughout this patent document, is defined broadly to mean inclusion 
without limitation. 

Exemplary control center 105 illustrates a central area that is commonly operator 
manned (not shown) for centrally monitoring and for centrally controlling the three 
exemplary process stages. A first process stage includes three raw material grinders 
1 10a to 1 10c that operate to receive a "feed" of raw material core and to grind the same, 
such as using a pulverizer or grinding wheel, into smaller particles of raw material. The 
term "or," as it is used throughout this patent document, is inclusive, meaning and/or. 
The second process stage includes a washer HOd that operates to receive the ground 
raw materials and clean the same to remove residue from the first stage. The third 
process stage includes a pair of separators 1 lOe and 1 lOf that operate to receive the 
ground and washed raw materials and separate the same, such as into desired minerals 
and any remaining raw materials. As this process facility is provided for illustrative 
purposes only and the principles of such are known, further discussion of the same~i7 
beyond the scope of this patent document. 

Exemplary control system 115 illustratively includes a global controller 120 and 
six local controllers 125a to 125f, each of which is implemented in software and 
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executable by a suitable conventional computer system (e.g., standalone, network, etc.), 
such as any of Honeywell, Inc.'s AM K2LCN, AM K4LCN, AM HMPU, AxM or 
like systems. Those skilled in the art will understand that such controllers may be 
implemented in hardware, software, or firmware, or some suitable combination of the 
same; in general, the use of computing systems in control systems for process facilities 
is known. 

Global controller 120 is associated with each of local controllers 125, directly or 
indirectly, to allow communication of information between the same. The phrase 
"associated with" and derivatives thereof, as used throughout this patent document, may 
mean to include within, interconnect with, contain, be contained within, connect to or 
with, couple to or with, be communicable with, cooperate with, interleave, be a 
property of, be bound to or with, have, have a property of, or the like. 

Global controller 120 monitors measurable characteristics (e.g., status, 
temperature, utilization, efficiency, cost and other economic factors, etc.) of associated 
processes 110, either directly or indirectly (as shown, through local controllers 125 
associated with processes 1 10). Depending upon the implementation, such monitoring 
may be of an individual process, group of processes, the facility as a whole, or 
otherwise. Similarly, local controllers 125 monitor associated processes 1 10, and, more 
particularly, monitor certain characteristics of associated processes 110. 

Global controller 120 generates, in response to such monitoring efforts, control 
data that may be communicated via local controllers 125 to associated processes 110 to 
optimize process facility 100. The phrase "control data," as used herein, is defined as 
any numeric, qualitative or other value generated by global controller 120 to globally 
control (e.g., direct, manage, modify, recommend to, regulate, suggest to, supervise, 
cooperate, etc.) a particular process, a group of processes, a facility, a process stage, a 
group of process stages, a sequence of processes or process stages, or the like to 
optimize the facility. Local controllers 125 operate to varying degrees in accordance 
with the control data to control the associated processes, and, more particularly, to 
modify one or more processes and improve the monitored characteristics and the 
facility. 

According to an advantageous embodiment, the control data may be 
dynamically generated using a lookup table defined and populated in accordance with 
the principles hereof, and such control data generation is based, at least in part, upon a 
given facility's efficiency, production or economic cost, and, most preferably, all three. 
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The lookup table may be populated and maintained on-line (e.g., at global controller 
120, at local controller 125, distributed within control system 1 15, etc.), off-line (e.g., 
standalone computer, network computer, etc.), or through some suitable combination of 
the same; likewise, the lookup table may be static upon population, be dynamic, or be 
modifiable, at least in part. 

The global controller 120 and the local controllers 125 may suitably use one or 
more such lookup tables to control processes 1 10 to conserve processing resources and 
increase the overall speed of control system 1 15. Control system 1 15 achieves a high 
level of both global and local monitoring, and cooperative control of associated 
processes 110 among controllers 120 and 125, by allowing the local controllers 125 to 
vary their individual or respective compliance with the control data. Varying degrees of 
compliance by local controllers 125 may range between full compliance and 
noncompliance. The relationship between global controller 120 and various ones of 
local controllers 110 may be master-slave (full compliance), cooperative (varying 
compliance, e.g., using control data as a factor in controlling the associated processes), 
complete disregard (noncompliance), as well as anywhere along that range. 

Depending upon the implementation and needs of a given facility, the 
relationship between global controller 120 and specific local controllers 125 may be 
static ( i.e., always only one of compliance, cooperative, or noncompliance), dynamic 
(i.e., varying over time, such as within a range between compliance and noncompliance, 
some lesser range therebetween, or otherwise), or varying between the same. One or 
more specific processes 110, and facility 100 as a whole, may be dynamically and 
cooperatively controlled as a function of local and global optimization efforts, and such 
dynamic and cooperative control is independent of the relationship between global 
controller 120 and specific local controllers 125, as described above. 

Turning to FIGURE lb, illustrated is a more detailed block diagram of one of 
the exemplary local controllers 125 that is associated with one or group of associated 
processes 110. Local controller 125 uses a single loop model predictive control 
("SL-MPC") structure that uses an efficient matrix prediction form in accordance with 
the principles of the present invention, as well as an analytical control solution map to 
reduce utilization of processing resources relative to conventional MPC technology. ~* 
According to the illustrated embodiment, local controller 125 receives as inputs, 
control/optimization specifications 130 (e.g., bounds, ranges, tolerances, control points, 
etc.) and feedback data 135 (e.g., output of associated process 110). 
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Control/optimization specifications 130 may be received from any of a number of 
sources depending upon the associated process or group of associated processes 1 10, an 
associated process facility or any other factor. For example, any of control/optimization 
specifications 130 may be received from an operator of a control center for the 
associated process facility, retrieved from a database or data repository, received from 
another associated controller (e.g., one or more local controllers 125, global controller 
120, or a suitable combination thereof), etc. 

Control/optimization specifications 130 include two types of variables: (1) a first 
variable ("MV") that may be manipulated, such as flow, feed, air blower, etc; and (2) a 
second variable ("DV") that cannot be manipulated and is a disturbance variable, such 
as bum rate, fuel quality per unit, etc. Feedback data 135 is a third variable ("CV") that 
is responsive to MVs and DVs, and is an output of associated process 110, such as 
pressure, temperature, etc. A sub-variable ("PV") of Feedback data 135 is indicative of 
the iterative response of the associated process 110 to monitoring and control by the 
local controller 125. Many, if not all, of such MVs, DVs and CVs represent measurable 
characteristics of associated process 110 that may be suitably monitored by local 
controller 125. 

Local controller 125 includes a dynamic prediction task with state estimation 
150, a local linear program/quadratic program ("LP/QP") optimization task 155, a 
dynamic control/optimization augmented range control algorithm ("RCA") 160 and a 
lookup table 165. Exemplary dynamic prediction task 150 receives CVs and operates to 
generate an array of multiple predictions (or dynamic unforced predictions) and, at 5 tau 
(response time close to end), an unforced prediction for values associated with 
associated process 110. The CVs represent feedback data 135 (e.g., inputs, outputs, 
etc.) associated with process 105, and dynamic prediction task 150 operates to accesses 
lookup table 165 and selects one or more values from the range of possible values, such 
selection being responsive, at least in part, to the received feedback data 135. A 
preferred method of using data structures, such as lookup table 165, or functionally 
equivalent dedicated circuitry, to maintain a range of possible values for one or more 
measurable characteristics associated with a process is disclosed and described in 
United States Patent Application Serial No. (Attorney Docket No. 120 25207), entitled 
"Process Facility Control Systems Using an Efficient Prediction Form And 
Methods of Operating The Same" and filed concurrently herewith, the disclosure of 
which has previously been incorporated herein by reference for all purposes as if fully 
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set forth herein. 

Exemplary local LP/QP optimization task 155 receives optimization 
specifications 140a and, in response to the unforced prediction, operates to generate, at 
5 tau, optimal values associated with associated process 1 10. 

A preferred method of performing the foregoing task is disclosed and described 
in United States Patent No. 5,758,047, entitled "Method of Process Controller 
Optimization in a Multivar.able Predictive Controller," which is commonly 
owned by the assignee of this patent document and related invention, the disclosure of 
which has previously been incorporated herein by reference for all purposes as if fully 
set forth herein. Most preferably, optimization specifications 140a are associated 
directly or indirectly, with an economic value of the output of associated process 1 10 
According to an advantageous embodiment, the unforced prediction may suitably be 
represented as a single variable and the LP/QP optimization task may be a linear 
determination of a minimum value or a maximum value, or a quadratic determination of 
a desired value. Exemplary dynamic control/optimization augmented RCA 160 
receives control specifications 140b and, in response to receiving the array of multiple 
predictions (from dynamic prediction task 150) and the optimal values (from local 
LP/QP optimization task 155), operates to generate control values, the MVs, that are 
input to associated process 1 10. An important aspect of exemplary local controller 125 
is the use of control/optimization specifications 140 and feedback data 135 to locally 
unify economic/operational optimization with MPC dynamically for a specific process 
or group of processes. 

Note the distinction between the foregoing discussion which introduces a very 
powerful multi-loop MPC embodiment having a well defined and dynamic 
interaction/interleaving relation among global and local controllers and the single loop 
controller embodiment described in United States Patent Application Serial No 
(Attorney Docket No. 120 25207), the disclosure of which has previously been 
incorporated herein by reference for all purposes. Those skilled in the art will 
understand the relationship among these embodiments and the applicability of the 
principles of the present invention. 

Turning now to FIGURE 2, illustrated is a flow diagram of an exemplar^ 
method (generally designated 200) for populating a data structure 165, shown as a 
lookup table, in accordance with the principles of the present invention (this discussion 
of FIGURE 2 makes concurrent reference to FIGURES la and lb). The phrase "data 
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structure," as the same is used herein, is defined broadly as any syntactic structure of 
expressions, data or other values or indicia, including both logical and physical 
structures. A data structure may therefore be any array (i.e., any arrangement of objects 
into one or more dimensions, e.g., a matrix, a table, etc.), or other like grouping, 
organization, or categorization of objects in accordance herewith. 

For purposes of illustration, a processor 205 and a memory 210 are introduced. 
Exemplary memory 210 is operative to store, or to maintain, lookup table 165, along 
with the various tasks/ instructions (generally designated 215) comprising method 200. 
Exemplary processor 205 is operative to select and execute tasks/ instructions 215 
which, in turn, cause processor 205 to perform the functions of method 200. 

To begin, processor 205 is directed through the execution of method 200 (e.g., 
manually (i.e., through interaction with an operator), automatically, or partially- 
automatically) to define a model 220 of at least a portion of at least one of the 
associated processes 110 (process step 225). Processor 205 is directed to store model 
220 in memory 210, preferably representing at least a portion of process 110 
mathematically. The mathematical representation defines one or more relationships 
among inputs and outputs of process 1 10. 

According to an advantageous embodiment, model 220 is defined using the 
following discrete state space model form: 

x k+l = Ax k + Su MD 

y k = Cx k + Du k(2) 

wherein x k , u k , and y k represent various states of modeled process 1 10, wherein A: is a 
time period and k+J is a next time period, and A, B, C, and D respectively represent 
measurable characteristics of modeled process 1 1 0 at any given time period. 

Processor 205 is directed to define a data structure, such as lookup table 165, 
having a plurality of accessible fields (process step 230). An exemplary source code 
embodiment for performing this definition is attached as Appendix A, and incorporated 
herein by reference as if fully set forth herein, and that is written in Pascal. Depending^ 
upon the needs of the particular implementation, the contents of such accessible fields 
may suitably be nulled, defaulted, or otherwise initialized or used. Memory 210, 
directed by processor 205, maintains lookup table 165, preferably representing, at least 
in part, an AB0I matrix 235 and a feedback vector 240. 
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According to an advantageous embodiment, AB0I matrix 235 and feedback 
vector 240 have the following respective definitions: 

\A B] 

0 I 



(3) 



(4) 



wherein / and 0 respectively and illustratively represent an identity matrix and an null 
matrix, for the purpose of this illustrative model, to maintain, or hold, MV constant 
(illustrated with respect to FIGURE 3). 

Processor 205 is directed to delineate mathematically a relationship among the 
above-given matrix 235 and vector 240 (process step 245), which according to an 
advantageous embodiment, has the following form: 

lA B\ 



'k+1 



0 J 



Z k = 



u 



(5) 



Processor 205 is directed to delineate mathematically a relationship among the above- 
given discrete state space model form and the Z vector 240 (process step 250), which, 
according to an advantageous embodiment, gives the following prediction form for any 
p interval, or point in the future: 



A r IA Bl P 

y(k +P )\k = [cd| q i . 



J k(6) 



Stated generally, use of Z vector 240 represents, or defines, mathematically, the 
relationship among the one or more inputs and outputs of modeled process 110. 

For a variety of purposes, as above-stated, for monitoring and for control of 
process 1 10, it is desirable to decrease utilization of processing resources. This may be 
accomplished, in part, through a recognition that certain characteristics of process 110 
are measurable (e.g., appraising, assessing, gauging, valuating, estimating, comparing^ 
computing, rating, grading, synchronizing, analyzing, etc.), whether or not such 
characteristics are dependent, independent, interdependent, or otherwise effected by 
other characteristics of the same process, a group of processes, a facility, a process 
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stage, a group of process stages, a sequence of processes or process stages, or the like. 
Many of these measurable characteristics have a range of possible values, which may or 
may not change, or vary, over time. It is desirable, in the present example, to determine 
an efficient prediction form ("EPF"), the range of values of which may suitably be 
maintained in lookup table 165. 

Processor 205 is directed to populate ones of the accessible fields 255 of lookup 
table 165 with a range of possible values of at least one measurable characteristic 
associated with at least process 110 (process step 260). An exemplary source code 
embodiment for performing this population is attached as APPENDIX B, and incorporated 
herein by reference as if fully set forth herein, and that is written in Pascal. According to 
the illustrative embodiment, it is desirable to have future predictions available, or 
precalculated, which may suitably be stored as an array of points within lookup table 
165. This collection of points may be referred to as PV-blocking, which may be given 
by the following form for any p, interval, or point in the future: 

yik+p^k 



A . 

Y(k+ pv—blocklng)y: 



P(k+p 2 )\k 

y(k+ Pm )\k 



(7) 



wherein i is the index for VV-blocking. The foregoing calculation may suitably be 
condensed into a product of EPF and Z*, which may be given by: 



y(k+pv-blocking)\k = [ EPF ] Z ^ (8 ) 

wherein EPF may be given by: 

~epf ± 

[EPF] = 



epf 2 
e P f m 



(9) 



wherein epfi is independent from the feed back information contained in the Z vector 
and may therefore be calculated in advance and given by: 
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A B 



•lP. 



epf . = \CD] 

1 \0 I 



(10) 



In short, exemplary processor 205 uses model 220 iteratively, or incrementally, 
to populate lookup table 165 with k possible values, thereby defining a range of values. 
A v* vector is formulated to conveniently calculate both Z k and yfi+pvbiockmgHk for 
different incremental k, which has the following form: 

Z , 





' AB ' 




v k+l = 


01 






EPF 





A 



Y(k+ pvbl ockin g)\k 



(11) 



Turning momentarily to FIGURE 3, illustrated is an exemplary two-dimensional 
graphical representation of MV and PV curves in accordance with a use of lookup table 
165 in accordance with the control system 100 of FIGURES la and lb and the 
principles of the present invention. It should be noted, that FIGURES la, lb, 2, and 3, 
along with the various embodiments used to describe the principles of the present 
invention in this patent document are illustrative only. To that end, alternate 
embodiments of model 220 may define any particular process, a group of processes, a 
facility, a process stage, a group of process stages, an interrelationship among, or a 
sequence of, processes or process stages, or some suitable portion or combination of any 
of the same. It should be further noted that a matrix structure was chosen for the EPF in 
this embodiment, however, alternate embodiments may use any appropriate data 
structure or dedicated circuitry to create a suitably arranged lookup array, or table, or 
the like. Such data structures and dedicated circuitry may be populated off-line, on-line 
or through some suitable combination of the same; likewise, such populated data 
structures and dedicated circuitry may be static, dynamic, modifiable, centralized, 
distributed, or any suitable combination of the same. 

Those of ordinary skill in the art should recognize that the computer system 105 
described using processor 205 and memory 210 may be any suitably arranged hand- 
held, laptop/notebook, mini, mainframe or super computer, as well as network 
combination of the same. In point of fact, alternate embodiments of computer system 
205 may include, or be replaced by, or combined with, any suitable circuitry, including 
programmable logic devices, such as programmable array logic ("PALs") and 
programmable logic arrays ("PLAs"), digital signal processors ("DSPs"), field 
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programmable gate arrays ("FPGAs"), application specific integrated circuits 
("ASICs"), very large scale integrated circuits ("VLSIs") or the like, to form the 
processing systems described and claimed herein. To that end, while the disclosed 
embodiments require processor 205 to access and to execute a stored task/instructions 
from memory to perform the various functions described hereabove, alternate 
embodiments may certainly be implemented entirely or partially in hardware. 
Conventional processing system architecture is more fully discussed in Computer 
Organization and Architecture, by William Stallings, MacMillan Publishing Co. (3rd ed. 
1993); conventional processing system network design is more fiilly discussed in Data 
Network Design, by Darren L. Spohn, McGraw-Hill, Inc. (1993); and conventional data 
communications is more fully discussed in Data Communications Principles, bv R. D. 
Gitlin, J. F. Hayes and S. B. Weinstein, Plenum Press (1992) and in The Irwin 
Handboo k of Telecommunications , by James Harry Green, Irwin Professional 
Publishing (2nd ed. 1992). Each of the foregoing publications is incorporated herein by 
reference for all purposes. 

APPENDIX A 

{ System name: SPID } 

{ BEGIN HDR } 

{ FILE: /lnk/spid/prv/src/sp_data.i.s } 

{VER/REV: 01,01 } 

{DATE: 05/14/98 } 

{ Honeywell Trade Secret - Treat as Honeywell PROPRIETARY } 

{END HDR } 

«*:PCHP:} 
unit SP DATA; 

interface 

uses mc_data, urv_data; 
{:PCHP:*)} 

CONST 

SPIDVERSION - 100.0000; {current version} 
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i iidA__il OIK 


= 10; 


maxnden 


= 5; 


rnax_nurn pin 


= 6; 


maxdelaydef 


= 200; 


URV tol 


= 0.001; 


UN 




Urr 




WARM 




FORCE 




AUTO 





= l; 

= 0; 

= 11; 
= 1; 



= 0; 



D_GRT_THAN_MAXD = 20 1 • 

BAD_NUM_OR_DEN = 20 2- 

BAD DEN INTEGRATOR =203; 

NEGATTVE_DEN_COEF =204; 

ZERO_GAIN =205 . 

ALLOC_S_Z_ERR = 20& 

BADMAXDINT = 207; 

BAD _T_CON = 2Q8 . 



TYPE 



sp_block_arr = array [1 ..max_n_blk] of single; 
sp_block_arr_ptr = A sp_block_arr; 

coef_arr = array [1 ..max_n_den] of single; { Predictor model F-polynom 

coefficients } 

coefarr_ptr = A coef_arr; 

coef_arrs = array [l..max_n_den * max_num_pin] of single; { Predictor 
model F-polynom. coefficients } 

sp_pin_arr_i = array [1 ..max_num_pin] of integer; 
sp_pin_arr_s = array [1 ..max_num_pin] of single; 

sp_pool_data_ptr = ^pjoo^data; 
sp_cds_data_ptr = A sp_cds_data; 
spcdsdata = record 
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{ 

- Section 1: Model Data Input Seciton (must have init value/BLD_visible) 
„ - } 



{ No. of DVs } 



{continuous model} 
{spid_cds:CAL} 
n_dv : single; 

{used as integer} 

{spid_cds:CALCl(1..6)} 

nnum : sp_pin_arr_s; { No. of B-poly coefficients } 

{used as integer} 
{spid_cds:CALC2(l ..6)} 

n_den : sp_pin_arr_s; { No. of F-poly coefficients } 

{used as integer} 

{spid_cds:CONVl(1..30)} 

num : coef_arrs; { Predictor model 

B-polynom. coefficients } 

{spid_cds:CONV2(l ..30)} 

den :coef_arrs; { Predictor 

F-polynom. coefficients } 

{spid_cds:CD(1..6)} 

delay : spjpin_an\s; { Dead-time (continuous) } 

{used as integer} 
{spid_cds:CC(1..6)} 

q s : S p_pin_arr_s; { Steady-state gains, used by 

controller } 

{spid_cds:CMD(1..6)} 

max_delay : sp_pin_arr_s; { Max dead-time (continuous) } 

{used as integer} 
{spid_cds:CPROF} 

perf rat : single; { Performance ratio - feedback } 
{spid_cds:CTREND3} 

Three Tau : single; { Three Tau 



model 
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(continuous) } 

{used as integer} 

{discrete model} 
{spid_cds:DEFl(1..6)} 

0 :s P-Pin_arr_s;{#ofnumcoef (discrete) 

} {used as integer} ' 

(spid_cds:DEF2(1..6)} 

(discrete* I Y a ■ '' SPJ?injUT - s > < # of dnum coef 

discrete) } {used as integer} 

{spid_cds:DEMNDINl (1 ..30)} 
b_z 

(discrete) } -coefaiTs; { num coefs 

{spid_cds:DEMNDIN2(l ..30)} 
f_z 

(discrete) } : coef_airs; { dnum coefs 

{spid_cds:DDD(1..6)} 

d : sp_pin_arr_s; { Dead-time (intervals) } 

{used as integer} 

{spid_cds:DDFPNT(l ..6)} 

max_d_int . S p_pin_arr_s; { maxd I T_con; (intervals)} 

{used as integer} «=rvais; ; 

{spid_cds:TSLOAD} 

T_con . sul gi e . 

• t , , . { Control execution 

interval (minutes) } 

{spid_cds:NUMSTIL} 



{- 



- Section 2: Variables used by Pascal/CL/LCN GUI 



limit.} 



{spid_cds:YY_L} 
y_L_ent 



: single; 



{Entered CV low 
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limit.} 



{spid_cds:YY_H} 

y_H_ent : single; 



{spid_cds:Y_L> 

y_L 



single; 



{Entered CV high 



{ CV low limit 



{spid_cds:Y_H} 

y_H 



single; 



{ CV high limit 



: single; 



value 



range } 



service } 



{spid_cds:YY} 

y 
} 

{spid_cds:Y^NOTOK> 
y_not_ok : single; 



{spid_cds:YVALl} 

first on : single; 



{spid_cds:YSEL} 

y hat : single; 

{spid_cds:Yl} 

yjast : single; 

{spid_cds:YVAL2} 

yFiltJBias : single; 

updating,Pool? } 

{spid_cds:YVAL3} 
Filt_Const : single; 

Pool? } 

{ spid_cds: YLDNEW} 
y Filt Ramp : single; 
{ spid_cds: YLDTGT } 
y_Filt_Ramp2 : single; 
{spid_cds:YIELD} 
Filt_opt_Const : single ; 



{ Current source 



{ 1: -> y value is uncertain/out of 



{ first time after pv 



{ Last good y value 



{ Filtered y-bias for bias 



{Bias correction Filter Factore, 
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{spid_cds:YIELDINT} 
Filt_ppt_Const2: single ; 
{spid__cds:Y2} 

y_last_BAD : single; 

BAD, 0=> Good ,P} 

{spid_cds:TYPE} 

y_type : single; 

LCN GUI uses it } 

{spid_cds:YJPRI} 

yJJFPO : single; 

{spid_cds:Y_NEWVAL( 1 ..2 1 )} 

yJJFP : arTay [0..20] of single; { 

prediction} 

{spid_cds:UNDIV01} 
u_L_ent : single; 



{ 1=> Last y is 



{ 0: stable, 1 rintegrator, 



limit.} 



limit.} 



{ spid_cds:UNDI V02 } 

u H ent : single; 

{spid_cds:U_L} 

u_L : single; 



Unforced 



{Entered MV low 



{Entered MV high 



{ MV low limit 



{spid__cds:U_H} 

uJH : single; 



{ MV high limit 



{spid_cds:LOWLIM} 
du_L : single; 



{ MV move low limit 



{spid_cds:HIGHLIM} 
du_H : single; 



{ MV move high limit 



{spid_cds:UCL(1..6)} 

u : sp_pin_arr_s; { Current actual value 
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{spid_cds:Ul(1..6)} 

uml : sp_pin_arr_s; { MV at previous interval 

} 

{spid_cds:U_0} 

u_0 : single; { Output from controller 

} 

{spid_cds:DIS_NAME} 

ss_display : single; { PC use only, use in simulation. Flag: 0 

means display ss values } 

{spid_cds:STP} 

saved_3t_pr : single; { saved Three_Tau ( for integrator ) or 

Perfrat ( for stable ) } 

{spid_cds:FLAGS> 

setup flag : single; { if = force, call setup; if = auto, don't call 

setup except necessory } 

n_sJL : single; { No. of intervals for J.L.'s 

mods} {used as integer} 

{spid_cds:BLOLVL( 1 .. 1 0)} 

y_blk : spblockarr; { CV blocking values 



} 



} 



{spid_cds:BLOMNMIN} 

y_num_blk : single; { total # of blocking intervals 



{spid_cds:BLOTIPC} 

y_num_start : single; { starting # of actual blocking 

intervals } 

{spid_cds:BLOMIN} 

y_num_end : single; { actual # of blocking intervals 

} 

{spid_cds:DIFPRXX(l ..5)} 

dx : coef_arr; { used in spidScalc} 

{spid_cds:ZOOMDEV(l ..5)} 

z : coef arr; { used in spidScalc} 
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{spid_cds:NUMSAMP} 
n_mvblk : single; 

time} {used as integer} 

{spid_cds:BLRQ( 1 .. 1 0)} 
mvblk 

{rmpccds.MODFACT} 
controller_mode : single; 
0-oflf; 1 -on; 1 1 -warm} 

{spid_cds:VVP} 

VerSWn ■ single; (save the release version in 



{ mv block sampling 

: sp_block_arr; { used in spidScalc} 
{only on DLL for now, 



cds} 



kt} 



{spid_cds:INITP ASS } 
initpass : singIe . 

{spid_cds:COUNTER} 
predi counter : single; 



{ Initialization pass counter 



{useful in debugging, the same as 



{ Current alarm 

: sp_block_arr; 

array [1 ..max_n_den * 



{spid_cds:Dl} 

now_alarm : single; 

indicators } 

{spid_cds:QANAPNT(l .. 1 0)} 
QMapcds 

{spid_cds:ENTMPO 1 1 (1 .. 1 50)} 
epflcds 

max_n_den * max_num_pin] of single; 

{spid_cds:ENTMPO 1 2( 1 ..300)} 

single- CPf2CdS ' ^ fL maX - n - blk * ma *-"_den * max_num_pin] of 

{spid_cds:DY(L.51)} 

. , ^ : array [0-.max_delay defl of 

single; { Output from controller } ~ — r 

{spid_cds:CDPNT} 

myl °° : s P-Cds_data_ptr;{ Initialization pass 

counter} K 
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(spid_cds:SPARE(1..30)} 

spares : array [1 ..30] of single; 



{ 

- Section 3: CL/LCN GUI Only 



10 



{spid_cds:TIMERQST} 
time_diff : single; 

CPU time in milliseconds } 



{ The amount of 



15 



20 



end; 

{KYS - save the whole QMap or QMapl which is the 1st row of QMap ? } 
sp _pool__data = record 

: sp_block_arr_ptr; 

: integer; 

: heap_array_s_ptr; 
: MatrixJType; 
: MatrixJType; 
: Matrix_Type; 



y_end 
n_s 
s_z 
epfl 
ep£2 
QMap 



end; 

{(*:PCHP:> 
25 implementation 
end. 

{:PCHP:*)} 



{ end sp_data } 
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APPENDIX Tt 

{ System name: SPID j 
{ BEGIN HDR j 

{ LCN FILE: /lnk/rmpc/prv/src/sp_proc.s } 
{ PC FILE: sp_proc.pas j 
{ DATE: 05/20/98 j 

{ Honeywell Trade Secret - Treat as Honeywell PROPRIETARY } 
{END HDR j 

{(*:PCHP:> 

unit SP_PROC; 

{$N+,E- 8087 IEEE floating-point, no emulation} 

{$A+,B-,G+ Word align, Boolean shortcut, protected mode} 
interface 

uses sp.data, mc_data {:PC_USE:}, urv_data {:}, mcjib, mc_rca, mem_mnag; 
procedure sp_ctfcheck(model_no 



i__no : integer; 


n_dv ; 


integer; 


var n_num 


: single; 


varnden 


single; 


var num 


: coef arrs; 


var den 


: coef arrs; 


var G_s 


: single; 


delay 


: single; 


var maxdelay 


: sp_pin_arr_s; 


var maxd _int 


: sp__pin_arr_s; 


T_con 


: single; 


var setup_flag 


: single; 


ThreeTau 


: single; 


perf rat : single; 


y_type 


: single; 


var saved_3t _pr 


: single; 
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var call_setup : single; 

var status : single); 

procedure getJJtau (n_den : integer; 

var den : coef_arrs; 

var ThreeJTau : single; 

var status : single); 

procedure sp_blocker ({ Inputs - } 

delay_int : single; 
{ Dead-times. } 

n_sJL : single; 
{ Number of intervals for J.L.'s mods. } 

perfrat : single; 

{ FF cntl/FB cntl response ratio.} 
{ Outputs ~ } 

var y_blk : spjblockjarr; 
{ CV block intervals (1 is current).} 

var yjriumblk : single; 
{ No. of block intervals per CV.} 

var Rtc : single; 

{ the beginning of control horizon} 

var status: single 
{ Status of blocker calcs}); 

Procedure spid_calc_S_z (n_B ♦ single; 

var B__z : coef_arrs; 

n_F : single; 

var F_z : coef_arrs; 

kp • integer; 

var S_z : heap_array_s_ptr; 

var status : single); 



procedure sp_ctrans (model_no : integer; 
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n_num : single; 

n_den : single; 

var num : coef_arrs; 

var den : coefarrs; 

G_s : single; 

T_con : single; 

{output} 

var n_b : single; 

var n_f : single; 

var b_z : coefarrs; 

var f_z : coef_arrs; 

var status : single); 



procedure sp_dt£2epfp (n_b 



: single; 
n_f : single; 

var b z : coef arrs; 

var fz : coefarrs; 

y_num_end : integer; 
var y_end : sp_block_arr; 

: single; 
: integer; 

: integer; 



delay_int 
lenunit 
norder 
{output} 
var epfl 
varepf2 
var status 



: Matrix_Type; 
: Matrix Type; 
: single); 



procedure new alloc urv (var urv : urv_set; 

totjnum_cv, {will include num of combined constraints for 
DQP} 

nummv : integer; 

callerlD : integer; {see URV JD ATA. PAS 

for definition} 
lenunit {:LCN_USE: : Sunitidentifier; :} 
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{:HPJJSE: : integer; :} 
{:VAXJJSE: : integer; :} 
{:GHSJJSE: : integer; :} 
{ :PCJJSE: } : integer; { : } {Points unit } 
single); 

(*:LOCAL:*) 



var status 



procedure new_release_urv(var urv : urv set; 



Sunit identifier; :} 



lcnunit {:LCNJJSE: 



{:HPJJSE: : integer; 



{:VAXUSE: : integer; :} 



{:GHSJJSE: : integer; 



{ :PC_USE: } : integer; { :} { Point's unit } 

var status : single); 



(*:LOCAL:*) 



procedure sp_h2solut (n_s 

var n_mvblk : single; 
var mvblk 
var s z 
y_num_end 
var yend 



intervals (1 is current). } 



delay int 
y^type 
G_s 
n sJL 
lcnunit 

overridetol: single; 



: integer; 

sp_block_arr; 

: heap_array_s_ptr; 

: integer; 

: sp_block_arr; { CV block 

: single; 
: single; 



: single; 



: single; 
integer; 
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var QMap : Matrix Type; 

var status : single 

); 

implementation 
{:PCHP:*)} 



1 0 {&&&8L&&S1SC&&&&&&&&SC&&&C&&&&&&&& 





nrocedure so ctfcheclcf model no 


. integer, 




n Hv 

11 U V 


. integer, 




VcU 11 11U11I 


. single, 


15 


var n den 


• ^mgie, 




var num 


. uuci__arrs. 




var den 


: coef arrs; 




var G_s 


: single; 




delay 


: single; 


20 


var maxdelay 


: sp_pin_arr_s; 




var max d_int 


: sp_pin_arr_s; 




T con 


: single; 




var setupflag 


: single; 




Three_Tau 


: single; 


25 


perf_rat 


: single; 




y_type 


: single; 




var saved_3t__pr 


: single; 




var call_setup 


: single; 




var status 


: single); 


30 


var 

i, j, totalmaxdint, first non zero ind: integer; 





begin 



if status = 0 then 
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begin {status = 0} 

{ if (callsetup = force) or the condition below is satisfied, call 

setup } 

if( call_setup = AUTO ) and 

( ( setup_flag = FORCE ) or 

( (y_tyP e = 0) and (perf_rat o 

saved_3t_pr) ) or 

( (y_type = 1) and (Three_Tau o 

saved_3t_pr) ) ) then 

call_setup := FORCE; 

if call_setup=FORCE then 

begin {callsetup} 



do 



)do 



while (nden > 0 ) and ( den[round(n_den)] = 0 ) 

n den := n_den - 1; 
while (n_num > 0 ) and ( num[round(n_num)] = 0 

njnum := n num - 1 ; 

if (delay > max_delay[modelno]) then 

status := D GRT THAN MAXD 

else if T_con <= 0 then 

status := BAD T CON 



else if (n_num > max n_den) or (n_den > 

max n den) or (n_num >= n_den) 

or (n_num < 1) or 

(n_den < 2) then 
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ranges} 



status := BAD NUM OR DEN; 

if status = 0 then 

begin {n_num and nden are inside of 



0 



n_den ) and 



firstnonzeroind + 1 ; 



first non zero ind := 1 ; 

while ( first_non__zero_ind <= 

( den[first_non_2ero_ind] = 0 ) do 
first non zero ind •= 



roimd(max_delay[model_no] / T_con); 



max_dJnt[model_no] 



total max d int + 



total_max__d_int := 0; 
forj := 1 to(n_dv+l)do 
total_max_d_int 



max delay def ) then 
BAD MAXD INT 



round(max_d_i 
ntQ]); 

lf ( total_max_d_int > 
status 



)or 



BAD_DEN_INTEGRATOR 



else if ( first_non_zero_ind > n den 

( first_non_zero_ind > 3 ) then 
status 
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(first_non_zero_ind+l) to round(n_den) do 



den[first_nonjzero_ind] ; 



den[first_non_zero_ind] ; 



den[first_non_zero_ind] := 1; 



round(n_den) do 



then 



:= NEGATIVEDENCOEF; 



ZERO GAIN 



round(n__num) do 



num[j] := num[j] / num[l]; 



num[l]; 



else 



begin 
for 



den[j] := den[j] / 



G_s := G_s 



for i := 2 to 
if den[i] < 0 
status 



end; 

if status = 0 then 

if num[l] = 0then 
status 



else 



begin 



for j := 2 to 



G_s := G_s / 
num[l] := 1; 



end; 
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ranges} 
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end- {n_num and n_den are inside of 



end; {call_setup} 
end; {status = 0} 

end; 



procedure get_3tau ( 



n den 

den 

ThreeJTau 

status 

var 

1 : integer; 

mt : single; 

begin 

if status = 0.0 then 

begin {status = 0.0} 



: integer; 

var 



: coef arrs: 



: single; 
single); 



var 



var 



Three Tau := 0; 
fori := 1 ton_den-l do 
begin {for} 

ifi = lthen 
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mt := den[2] 
else ifi = 2then 

mt := sqrt(den[3]) 
else ifi = 3 then 

mt:=exp(ln(den[4])/3) 

else 

mt := sqrt( sqrt( den[5] ) ); 

if mt > Three JTau then 
Three Tau := mt; 



end; {for} 



Three Tau := Three_Tau * 3; 



end; {status = 0.0} 



end; 



&&&<! 

Procedure spid_calc_S_z( n_B : single; 



var B_z : coef_arrs; 
n_F : single; 
var F_z : coef_arrs; 
kp : integer; 

var S_z : heap_array_s jptr; 

var status : single 
);(* -.LOCAL:*) 



Var 



loc_y : array[-n_F_max_dim..01 of single; 

total_B, term 1 , term2 : single; 
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i> i_kp : integer; 

Begin {spid_calc_S_z} 

if status = 0.0 then 

begin {status = 0} 

{Calc total_B to save calculations later.} 

total B := 0.0; 

for i := 1 to round(nB) do 

total B := total B + B_z[i]; 

{Initialize Ioc_y array to 0} 
for i := 0 to round(n_F) do 
loc_y[-i] := 0.0; 

for i_kp := 1 to kp do 
begin {ijcp} 

{First, shift ("age") loc_jy array} 
for i := - round(n_F) to -1 do 
loc_y[i] := loc_y[i+l]; 

{Then, find terml - the contribution of B_z} 
if ijkp > round(n B) then 
terml :=total_B 

else 

begin 

terml :=B_z[l]; 
for i := 2 to ijkp do 

terml := terml + B_z[i]; 

end; 



{Next, calculate term2 - the contribution of F_z} 
term2 := 0.0; 

for i := 1 to round(n_F) do 
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term2 := term2 - F_z[i] * loc_y[-i]; 



{Finally, put current y into loc_y[0]} 
loc_y[0] :=terml +term2; 
s_z A [i_kp] := loc_y[0]; 



end; {i_kp} 
end; {status = 0} 
End; {spid _calc_S_z} 



{&&&&&&&&&^ 

procedure sp blocker ( 
{ Inputs ~ } 

delay_int : single ; { Dead-times. 

} 

n_sJL : single ; { Number of intervals for J.L.'s 

mods. } 

perf rat : single ; { FF cntl/FB cntl response 

ratio. } 

{ Outputs -- } 

var y_blk : sp_block_arr ; { CV block intervals (1 is current). } 

var y_num_blk : single ; { No. of block intervals per CV . 



} 



var 



var Rtc: single ; { the beginning of control horizon } 
var status: single { Status of blocker cedes } 

); 



Rte, Rti : single; 

i, nkvi : integer; 

Iv temp : cv block arr; 
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: cv block arr; 



begin { Start blocking calculations } 
if status = 0.0 then 

begin {status = 0} 

{for SISO, Rti = Resp time = n_sJL } 
Rti := n_sJL + delayint; 
Rtc := minr(l, perfrat) * njJL + delayint; 
Rte:=(Rti + Rtc)/2; 

GetPVblk( Kvjemp, 

Ivtemp, 
nkvi, 

round(delay_int), 

Rtc, Rte, 

10, 

1, 

1, 

false, 
1.0); 

for i := 1 to nkvi do 

y_blk[i] := Kv_temp[i]; 

ynumblk := nkvi; 
end; {status = 0} 
end; { End blocking calculations } 



&&&&&&&&&<&&&&&&&&&& 

procedure sp_ctrans (model_no . i nte ger 

n num 

single; 
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single; 



coefarrs; 



coef arrs; 



single; 



nden : 

var num : 

var den : 

G_s : single; 

T_con : single; 

{output} 

var n_b : 

var n_f : single; 

var b_z : coefarrs; 

var f_z : coefarrs; 

var status : single); 



var 

nnumi, n deni : integer; 

sqtau2, taulmul4, tmidl, tmid2, tmid3, tmid4, tmid5, tmid6, tmid7, tmid8, a, b 
: single; 

numptr, denptr, b z_ptr, f_zjptr : coef_arr _ptr; 

begin { ctrans } 

if status = 0.0 then 

begin {status = 0} 

nnumi := round(nnum); 

n_deni := round(nden); 

if (njiumi = 1) and (n deni = 2) and (den[l] = 0) then 
begin { 1/S } 

n_b := 1 ; 
n_f:=l; 

b_z[l] :=G_s*T_con; 
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else if (nnumi = 1) and (n_deni = 2) and (denfl] = 1 ) then 
begin {I/(tl*S + l)> 

n_b:= 1; 
n_f~l- 

b_z[l] :=G_s * (1 -exp(-T_con/den[2])); 
Lz[l] := - exp(-T_con / den[2]); 
end {l/(tl*S+l)} 

else if (n_numi = 1) and (n_deni = 3) and (denfl] = 1) then 
begin { l/( tl * S^2 + t2 * S + 1) } 

n_b:=2; 
n_f:=2; 

sqtau2 := sqr ( den[2] ); 
taulmul4 :=den[3] * 4; 
tmid 1 .= sqrt(abs( sqtau2 - taul mu!4 )); 

if(tmidl <eps)then 
begin 

tmid2 := T_con / sqrt(den[3]); 
tmid3 := exp (-tmid2); 
tmid2)); b - z n]:=G_s*(l- tal id3*(l + 

tmid2 - 1); b - Z[2J ' = °- S * toid3 * + 

f_z[lj := (-2) * tmid3; 
f_z[2] := tmid3 * tmid3; 

end 

{kys - think about sqtau2 = taulmuI4} 
else if ( S qtau2 < tau 1 mu!4 ) then — 
begin {Tau2*2 less then (4* Taul ), w0 



tmid8 := den[3] * 2; 
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tmidS + tmid3 ) ); 



tmid4 * ( tmid3 - tmid5) ); 



Taul ), wO } 



Taul ), ab } 
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tmid2 := tmidl / tmid8; 

tmid7 := tmid2 * T_con; 

tmid5 := cos(tmid7); 

tmid6 := sin(tmid7); 

tmid3 := den[2] / tmidl * tmid6; 

tmid4 := exp ( - den[2] / tmid8 * 

b_z[l] := G_s * ( 1 - tmid4 * ( 

b_z[2] := G_s * ( sqr(tmid4) + 

f_z[l] := -2 * tmid4 * traid5; 
f_z[2] := sqr( tmid4 ); 
end { Tau2 A 2 less then ( 4 * 



else if ( sqtau2 > taulmul4 ) then 

begin { Tau2 A 2 greater then ( 4 * 

a := (den[2] - tmidl) * 0.5 / den[3]; 
b := (den[2] + tmidl) * 0.5 / den[3]; 
tmid2 := tmidl / den[3]; {b-a} 
tmid3 := exp ( - a -* T_con ) ; 
tmid4 := exp ( - b * T_con ) ; 
b_z[l] := G_s * (b * (l-tmid3) - a * 



(l-tmid4))/tmid2; 

tmid3 - b * (l-tmid3) * tmid4) / tmid2; 

T_con ); 
Taul ), ab } 



end; 



b_z[2] := G_s * (a * (l-tmid4) : 

f_z[l] := - tmid3 - tmid4; 

f_z[2] := exp ( - den[2] / den[3] 

{ Tau2 A 2 greater then ( 4 
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end { l/(tl *S A 2 + t2*S+l)} 

else 

begin 

{(*:PC_ONLY:} 

writelnCCtrans cannot do it now'); 

{:PCJ3NLY:*)} 

{status := 100; } 

end; 

end; {status = 0} 

end; { ctrans } 

{&&&&&&&&&&&&&&&&&&&& 
&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 

{ [epfl , epG]=dtf2epfp(mn,md,pv_end-mdt,0) } 
procedure sp_dt£2epfp ( 

nb 

: single; 

n_f 

: single; 

var 

b_z : coef_arrs; 

var 

f_z : coef_arrs; 

y_num_end : integer; 

var 

y_end : sp_block_arr; 



delay_int 



: single; 
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lcnunit 



: integer; 



n order 



integer; 



{output} 

epfl 
ep£2 
status 
var 

pvrow 

nl, get amount, i, j, il, jl 

code : single; 

gen_ptr : r anyptr; 

begin { sp_dt£2epfp } 

if status = 0.0 then 

begin {status = 0} 
code := 0.0; 



; Matrix_Type; 
MatrixJType; 



single); 



Row: 



integer; 



var 



var 



var 



for i ;= 1 to n_order do 
begin 

epfl.mat[i] A [l] :=-f_z[i]; 
for j := 2 to n_order do 
if 0=1+1) then 

epfl.mat[i] A (j] := 1 

else 

epfl.mat[i] A [j] - 0; 

end; 
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if status = 0.0 then 

begin {no error after sp_tf2obsv} 

for i := 1 to n_order do 

epfl.mat[i] A [epfl.m] :=b_z[i]; 

for j := 1 to n_prderdo 

epfl.mat[epfl.n] A G] := 0.0; 

epfl.mat[epfl.n] A [epfl.m] := 1; 

forj := 1 toepfl.ndo 

pvrow[j] := epfl.mat[l] A [j]; 

i:=l; 

for j := 1 to round(y_end[y_num _end]-delay_int) 
begin 

Pre VectMult (pvrow, epfl , 

for j 1. := 1 to epfl .m do 

pvrowjjl] := 

if ( j = round(y_end[i]-delay_int) ) 

i :=i + 1; 

end; 

epfl.n := epfl.n -1; {taking out the last row, but 
leaving the memory } _^ 

end; {no error after sp_tf2obsv} 
end; {status = 0} 
end; { sp__dtf2epfpcal } 



ep£2.mat[i] A ); 



epf2.mat[i] A 01]; 
then 
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— } 

{Procedure to allocate space for U,R, and V. Called from setup } 
{and from rca_body (in order to get second URV for PreRCA) } 

procedure new alloc urv ( varurv : urv set; 

tot_num_cv, {will include num of combined constraints for DQP} 
num_mv : integer; 

callerlD : integer; {see URVJDATA.PAS for definition} 

lcnunit {:LCN_USE: : $unit_identifier; :} 

{:HP_USE: : integer; .} 
{:VAX_USE: : integer; .} 
{:GHS_USE: : integer; .} 
{:PC_USE:} : integer; {:} { Point's unit } 

status : single); (*:LOCAL:*) 



var 



var 



n,m,i : integer; 
code : single; 

begin {new alloc urv} 
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code := 0.0; 
with urv do 

begin {with urv} 

{Based on callerlD, calculate n and m.} 
if (callerlD = CID_RMPCT) or (callerlD 
CID_RMPCT_PRERCA) then 

begin 

{later, might determine different (smaller) size for 

PreRCA} 

n := maxi(tot_num_cv * max_cv_blk + num_mv * 
(max_mv_blk-l), tot num cv * 2 + nummv + 1); 

{for dynamic CV 
rows} {for dynamic CX rows} {for SS CV row + extra LP/QProws} 

m := num_mv * max mv blk; 

end 

else if (callerlD = CIDDQP) then 
begin 

n := tot num cv * 2 + num_m v + 1 ; 

{for SS CV row + extra 

LP/QProws} 

m := num_mv; 

end 

else if (callerlD = CID_SPID) then 
begin 

n := tot num cv; 
m := nummv; 

end 

else 

begin 

n := 0; 

m := 0; ^ 

end; 



n := mini(max_row_size, n); 
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{double check MC_Control/OptSS later, what if n > 
maxjro w_size??? } 

{store size allocated for later use in release_urv} 
n__alloc := n; 
m alloc := m; 

{grab space for U} 
for i := 1 to n do 

U.mat[i] := (*:loophole:*) Row_ptr ( new_mc_getmem(n 



* sizeof(URV_sd), 



* sizeof(URV_sd), 



lcnunit, code, status)); 

{grab space for R} 
for i := 1 to n do 

R.mat[i] := (*:loophole:*) Row_ptr ( newjmc _getmem(m 

lcnunit, code, status)); 

{grab space for V} 
for i := 1 to m do 

V.mat[i] := (* -.loophole:*) Rowjrtr ( new_mcjgetmem(m 



* sizeof(URV_sd), 



lcnunit, code, status)); 



U.n:=0; U.m:=0; 
R.n:=0; R.m:=0; 
V.n:=0; V.m:=0; 
R.mt := UpperTri; 



end; {withurv} 
end; { new al locjurv } 



{&&&&&&&&&&&&&&&&&&&&&& 
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procedure new_release_urv(var urv : urv set; 

lcnunit {:LCN_USE: : $unit_identifier; 

:} 

{rHPJJSE: : integer; : } 

{:VAX_USE: : integer; :} 

{:GHS_USE: : integer; :} 

{:PCJJSE:} : integer; {:} { Points unit } 
var status : single); (*:LOCAL:*) 

var 

n,m,i, getamount : integer; 
gen_ptr : ranyptr; 

begin {new_release_urv} 
with urv do 

begin {with urv} 

n := n alloc; 
m := malloc; 
{grab space for U} 
get_amount := n * sizeof(URV_sd); 
for i := 1 to n do 
begin 

gen__ptr := (*:loophole:*)r_anyptr(U.mat[i]) ; 
new_mc_freemem(gen_ptr, get amount, lcnunit, 



status); 



end; 



{grab space for R} 

get amount :=m * sizeof(URV_sd); 
for i := 1 to n do 
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begin 



gen_ptr := (* .loophole: *)r_anyptr(R.mat[i]); 
new_mc_freemem(gen_jptr, getamount, lcnunit, 



end; 

{grab space for V} 
get amount := m * sizeof(URV __sd); 
for i := 1 to m do 
begin 



gen_ptr := (*:loophole:*)r_anyptr(V.mat[i]); 
new_mc_freemem(gen_ptr, get_amount, lcnunit, 



status); 



status); 



end; 



end; {withurv} 
end; {new release urv } 




var 



nblk 



: integer; 



lastblk 



: integer; 



dense 



: integer; 



var 



mv blk 



: sp_block_arr 



); 



var 
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i : integer; 

begin {getmvblk} 

if nblk >= max_n_blk then 
nblk := max n blk; 



if (lastblk <= nblk) then 

begin { lastblk <= nblk } 

for i := 1 to lastblk do 
mv_blk[i] := i; 
nblk := lastblk; 
end { lastblk <= nblk } 

else 

begin { lastblk > nblk } 
mv_blk[l] := 1; 
for i := 2 to nblk do 

mv_blkp] := 0; 



1)* 



1] + 



for i — 2 to nblk do 
begin 

mv_blk[i] ;= mv_blk[i] + ( exp (dense * i / nblk) - 

( lastblk / ( exp(dense) - 1 ) ); 
if mv_blk[i] > ( mv_blk[i-l] + 1 ) then 
mv_blk[i] := round( mv_blk[i] ) 

else 

begin 

if i < nblk then 

mv_blk[i + 1] := mv_blk[i + 

( mv_blk[i - 1] + 1~T 
mv^blk[i])/2; 
mv_blk[i] := mv_blk[i - 1] + 1; 

end; 



WO 00/41045 PCT/US99/27726 

-50- 

end; 

end; { lastblk > nblk } 
end; {getmvblk} 

5 



10 



15 



20 



25 



{&&&&&&&&&&&&&&&&&&&&& 

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&} 
procedure sp_h2solut ( 

n_s : integer; 

var njnvblk : single; 



intervals (I is current). } 



30 



var mv blk : 
var s z 
y_num_end 
var y_end 

delay_int 

y_type 

G_s 

n_sJL 

lcnunit 

override tol: single; 

var QMap 

var status : single 

); . 



sp_block _arr; 

: heap_array_s_ptr; 

: integer; 

: sp_block_arr; { CV block 



: single; 
: single; 



: single; 



: single; 
: integer; 

: Matrix_Type; 



var 



return 



: single; 
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i, j , dense, n_mvblki, last_blk : integer; 



urv 



: urvset; 

built row 



Row; 




{Finds a value in the A matrix, which 



is never explicitly built as a } 



{complete matrix. The only parameters to A are i and j, the row and } 
{column indices in A. A does, however, access several of the parameters 

{that are passed to the main RCA procedure that calls A These "global" 



{parameters are 
{ 



{ pv_blk, 
column in our A,} 

{ 
{ 
{ 



> 



single array, index is row or 



value is index to entire matrix "S?" of step } 
response data j 

} 



{ delay_int - integer, value is dead time for each block } 



{ 

{ S_z 



{ n_S 

} 

{ 



} 

-- array of step response polynomials} 

} 

~ integer, value is number of values in 
step response polynomial \ 
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10 



15 



20 



25 



30 



Function A(i j : integer) : single; (*:LOCAL:*) 
var 

value : URV_sd; 

ind_S : integer; {index within step response polynomial} 

Begin {A} 

value := 0.0; 

{use the 5 tau value if process is unstable} 
{Q2J - verify} 

ind_S := round(y_end[i] - delay _int - j + 1); 
ifind_S>0then 
begin 

if ind_S <= n_S then 

value := S_z A [ind_S] 
else if ( y__type = 1 ) then {not stable} 

value := S_z A [n_S] + G_s * (ind_S - n_S) 

else 

value := S_z A [n_S]; 

end; 

{ 

if((i-j+l)>0)then 

value := S^z A [i-j+l]; 

} 

A := value; 



End; {A} 
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{start of sp_h2solut procedure body} 



begin {sp_h2solut} 

if status = 0.0 then 

begin (status = 0} 



dense := 3; 

n rnvblki := round(njnvblk); 

last_blk := round( minr(n_sJL, y_end[y_num_end] - delay Jnt) ); 

getmvblk (n_mvblki, lastblk, dense, mv^blk); 

new_alloc_urv ( urv, y_num_end, n_mvblki, CID_SPID, lenunit, 

status); 



if status = 0.0 then 
with urv do 

begin {with urv} 

return := successful; 
k := 0; 

{QMap = pinv(Hss) = V * (R \ LP) } 
{URVAddRow returns U, R, and V, no 

transpose } 



for i := 1 to y_num_end do 
begin 

for j := 1 to nmvblki do 

built_rowD] := A(i, 

round(mv_blk[j]) ); 

{ built_row[j] : = 

A_blk(ij);} 

return : = 

URVAddRow(U,R,V,k, built_row, n_mvblki); 
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U.mat[i] A D]; 

is U row into built_row, 

(R, built_row, k); 
{ 

built_row, QMap.mat[i] A ); 

transpose here 

the whole QMap on Am, 

} 

QMap.mat[l] A [i] + 
{ 

QMap.mat[l] A [i] / ( n_s/n_mvblk ); 



end; 

for i := 1 to R.n do 
begin 

for j : 

{R\U f , pass IT collom which 

builtrow is a used as collom} 

return ForwardBackSub 

PostVectMult (V, 

QMap is acturally QMap's 

We do not need to calculate 

but for offline, probably. 

QMap.mat[l] A [i] := 0; 
for j := 1 to k do 
QMap.mat[l] A [i] , := 

V.mat[l] A [j] * built_row[j]; 
QMap.mat[l] A [i] := 

} 

end; 

newrelease_urv(urv, lcnunit, status); 



= 1 to k do 
built_row[j] 



end; {with urv} 
njmvblk := n_mvblki; 
end; {status = 0.0} 
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end; {sp_h2solut} 



{(*:PCHP:> 
END. {of unit} 
{:PCHP:*)> 
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What is Claimed is : 

1- A computer system for use with a process facility having a plurality of 
associated processes, comprising: 

circuitry that is capable of maintaining a data structure having a plurality 
of accessible fields; and 

a processor, associated with said circuitry, that is capable of populating 
ones of said plurality of accessible fields of said data structure with a range of possible 
values of at least one measurable characteristic associated with at least one process of 
said plurality of associated processes. 

2. The computer system set forth in Claim 1 wherein said circuitry is 
capable of storing a task that directs said processor to populate said ones of said 
plurality of accessible fields of said data structure with said range of possible values. 

3. The computer system set forth in Claim 1 wherein said circuitry is 
further capable of maintaining a model of at least a portion of said plurality of 
associated processes. 

4. The computer system set forth in Claim 3 wherein said model includes a 
mathematical representation of at least a portion of said at least one process of said 
plurality of associated processes, said mathematical representation defining 
relationships among inputs and outputs of said at least one process of said associated 
processes. 

5. The computer system set forth in Claim 3 wherein said processor is 
capable of using said model iteratively to populate ones of said plurality of accessible 
fields of said data structure with said range of possible values of said at least one 
measurable characteristic. 

6. The computer system set forth in Claim 5 wherein said model includes at 
least one feedback variable representing, at least in part, an output of said at least one* 
process of said associated processes. 



7. The computer system set forth in Claim 6 wherein said processor 
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in 



populates at least one of said plurality of accessible fields of said data structure _ 
response to said at least one feedback variable of said at least one process of said 
associated processes. 



8. The computer system set forth in Claim 3 wherein said model includes 
manipulable variable. 



9. The computer system set forth in Claim 8 wherein said processor is 
capable of at least substantially maintaining a value of said manipulable variable during 
at least a portion of said iterative population of said ones of said plurality of accessible 
fields of said data structure. 

10. The computer system set forth in Claim 1, wherein said circuitry 
mamtains statically said range of possible values of said at least one measurable 
characteristic associated with at least one process of said plurality of associated 
processes. 

11. The computer system set forth in Claim 1 wherein said processor is 
further capable of using said range of possible values of said at least one measurable 
characteristic to predict an unforced response associated with said at least one process. 

1 2. A method of operating a computer system that is for use with a process 
facility having a plurality of associated processes, said method of operation comprising 
the steps of: 

maintaining a data structure having a plurality of accessible fields in 
circuitry associated with said computer system; and 

populating ones of said plurality of accessible fields of said data structure 
using a processor, that is associated with said circuitry, with a range of possible values 
of at least one measurable characteristic associated with at least one process of said 
plurality of associated processes. 

13. The method of operation set forth in Claim 12 further comprising the 
step of storing a task in said circuitry that is capable of directing said processor to 
populate said ones of said plurality of accessible fields of said data structure with said 
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range of possible values. 

14. The method of operation set forth in Claim 12 further comprising the 
step of maintaining a model of at least a portion of said plurality of associated processes 
in said circuitry. 

15. The method of operation set forth in Claim 14 wherein said model 
includes a mathematical representation of at least a portion of said at least one process 
of said plurality of associated processes, said mathematical representation defining 
relationships among inputs and outputs of said at least one process of said associated 
processes, and said method further comprises the step of using said model iteratively by 
said processor to populate ones of said plurality of accessible fields of said data 
structure with said range of possible values of said at least one measurable 
characteristic. 

16. The method of operation set forth in Claim 15 wherein said model 
includes at least one feedback variable representing, at least in part, an output of said at 
least one process of said associated processes, and said method further comprises the 
step of using said processor, in response to said at least one feedback variable of said at 
least one process of said associated processes, to populate at least one of said plurality 
of accessible fields of said data structure. 

17. The method of operation set forth in Claim 14 wherein said model 
includes a manipulable variable, and said method further comprises the step of at least 
substantially maintaining a value of said manipulable variable during at least a portion 
of said iterative population of said ones of said plurality of accessible fields of said data 
structure. 

18. The method of operation set forth in Claim 12 wherein said circuitry 
maintains statically said range of possible values of said at least one measurable 
characteristic associated with at least one process of said plurality of associated 
processes. 

19. The method of operation set forth in Claim 12 further comprising the 
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step of predicting an unforced response associated with said at least one process using 
said processor and said range of possible values of said at least one measurable 
characteristic. 

5 20. A data structure for use with a computer system associated with a 

process facility having a plurality of associated processes, said data structure comprising 
a plurality of accessible fields, ones of said plurality of accessible fields maintaining a 
range of possible values of at least one measurable characteristic associated with at least 
one process of said plurality of associated processes. 

0 

21. The data structure set forth in Claim 20 wherein one of said plurality of 
accessible fields are capable of being selected by said computer system. 
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